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An  Efficient  and  Accurate  Method  for  Calculating  and  Representing 
Power  Density  in  the  Near-Zone  of  Microwave  Antennas 

Richard  L.  Lewis  & Allen  C.  Newell 


An  algorithm  is  presented  for  calculating  near-zone  and  Fresnel- 
region  fields  in  front  of  microwave  antennas  from  discrete  numerical 
values  of  the  radiated  plane-wave  spectrum  (complex  far-field  pat- 
tern). That  is,  the  near  fields  are  calculated  by  numerically  inte- 
grating the  plane-wave  spectrum  representation  of  the  field.  The 
crux  of  the  analysis  consists  of  handing  a numerical  instability 
which  arises  from  integrating  discrete  data.  A criterion  is  developed 
for  limiting  the  integration  domain  in  order  to  exclude  highly  oscil- 
latory regions  of  the  integrand.  In  turn,  this  leads  to  restricting 
the  applicable  output  range  over  which  the  field  can  be  computed. 

With  the  numerical  instability  problem  thus  resolved,  fast  Fourier 
transform  techniques  are  used  to  assure  efficient  numerical  inte- 
gration over  a large  (but  restricted)  output  range.  The  results  are 
conveniently  presented  as  relative  power-density  contours  in  planes 
formed  by  the  longitudinal  coordinate  axis  and  one  transverse  coor- 
dinate axis.  The  algorithm  is  capable  of  extremely  high  accuracy, 
which  is  demonstrated  by  comparing  predicted  and  measured  near  fields 
for  two  distinct  antennas,  along  with  a comparison  against  an  exact 
theoretical  model.  In  the  case  of  circularly  symmetric  excitation 
models  for  electrically  large  antenna  apertures,  the  predicted  rela- 
tive near-zone  power-density  contour  plots  turn  out  to  be  a function 
of  just  the  relative  aperture  distribution.  Nomographs  for  obtaining 
absolute  near-zone  power  densities  are  presented  for  a few  typical 
aperture-di stribution  functi ons. 

Key  words:  antennas;  antenna  apertures;  near  fields;  Fresnel 

region;  plane-wave  spectrum. 


I.  Introduction 

There  are  a number  of  situations  related  to  safety  and  interference  where 
it  is  desirable  to  have  a reliable  and  concise  method  of  estimating  the  power 
density  levels  radiated  by  microwave  antennas.  We  present  a new  technique  for 
predicting  near-zone  field  intensities  using  plane  wave  scattering  matrix 
theory  to  produce  accurate  calculations  along  with  a highly  efficient  method 
of  graphically  representing  the  results.  Using  this  technique,  comparisons  of 
predicted  and  measured  near  fields  were  carried  out  for  selected  antennas  and 
found  to  be  in  excellent  agreement. 

The  input  to  the  calculations  is  the  far-field  pattern  (plane-wave  spec- 
trum) of  the  antenna  obtained  from  either  measured  near-field  data,  design 
specifications,  or  theoretical  analysis.  From  this  input,  programs  based  on 
the  plane-wave  theory  were  developed  to  calculate  the  field  along  the  x and  y 
axes  on  a series  of  planes  parallel  to  the  aperture.  Then  the  power-density 
values  in  the  x-z  and  y-z  planes  are  obtained  and  graphically  represented  on 
two  contour  plots.  The  ordinate  of  these  graphs  covers  the  range  ±4D  in  units 
of  D and  the  abcissa  extends  from  just  in  front  of  the  antenna  out  to  D2/\  in 
units  of  D2/A,  where  D is  the  antenna  diameter  and  A is  the  wavelength.  This 
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choice  of  scaling  means  that  a single  graph  will  apply  without  approximation 
to  all  antennas  with  the  same  D/A.  ratio  and  relative  aperture  distribution. 
Different  input  power  levels  can  be  accounted  for  by  scaling. 

Theoretical  results  are  presented  in  section  II  of  this  report  ignoring 
small  structure  effects  (which  frequently  do  not  scale  with  D/A),  producing 
graphical  presentations  which  just  depend  on  the  relative  aperture  distribu- 
tion regardless  of  the  D/A  ratio.  The  computational  method  is  presented  in 
section  III,  which  concludes  with  specific  examples  using  measured  data  to 
demonstrate  the  accuracy  of  this  technique. 


II.  A Quick  Method  for  Estimating  Near-Zone  Field  Intensities 
Radiated  by  Circular  Aperture  Antennas 


There  are  many  instances  where  an  estimate  is  needed  of  the  near-zone 
field  intensity  radiated  by  a parabolic  dish.  One  possible  need  is  to  evalu- 
ate hazards  to  personnel,  explosive  detonators,  or  sensitive  electronics.  A 
goal  of  the  U.S.  Environmental  Protection  Agency  (EPA)  is  to  characterize 
field  intensity  levels  in  the  near  field  of  typical  aperture  or  reflector 
antennas.  On  behalf  of  the  EPA,  the  U.S.  National  Bureau  of  Standards  (NBS) 
has  developed  techniques  for  predicting  antenna  near  fields.  A few  of  these 
predictions  are  compared  against  near-field  measurements  for  specific  antennas 
and  found  to  be  in  excellent  agreement.  Here,  we  wish  to  characterize  near- 
zone field  intensity  levels  for  a specific  class  of  antennas  using  these  tech- 
niques. 

Our  near-field-prediction  procedure  is  numerical  integration  of  the  plane- 
wave-spectrum  integral  representation  for  the  electric  field.  The  plane-wave 
spectrum  data  used  by  this  procedure  can  be  generated  either  from  measured 
data  or  analytical  models.  In  this  study,  we  consider  relative  aperture  dis- 
tribution functions  which  are  amenable  to  analytical  integration  to  obtain 
theoretical  plane-wave  spectra. 

We  have  chosen  for  study  a circular  aperture  with  a rotational ly  symmetric 
tapered  aperture  illumination.  We  wish  to  plot  relative  near-zone  power  den- 
sities as  a function  of  (normalized)  y and  z coordinates,  where  z is  measured 
along  an  axis  perpendicular  to  the  antenna  aperture  while  y is  measured  in  a 
plane  transverse  to  the  z axis.  The  expected  antenna  pattern  is  a function  of 
^sin6,  where  D is  the  antenna  diameter,  A is  the  wavelength,  and  0 is  the 
angular  position  of  the  observation  point  with  respect  to  antenna  boresight. 
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Consequently,  upon  normalizing  boresight  distance  z with  respect  to  D2/A  and 

transverse  distance  y with  respect  to  D,  we  find  that  a single  set  of  curves 

becomes  sufficient  to  describe  the  relative  near-field  antenna  pattern  for  a 

given  D/A  ratio  and  a given  relative  aperture  distribution.  Moreover,  if  the 

D2 

maximum  transverse  distance  y is  small  compared  to  t—  (we  assume  D > 30A 

■'max  A ~ 

and  y = 4D),  then  the  near-field  patterns  become  invariant  with  respect  to 
-’max  ’ 

D/A.  When  D < 30A  a separate  near-field  pattern  for  each  D/A  value  is 
required.  Fortunately,  the  patterns  change  slowly  so  one  can  interpolate 
between  just  a few  patterns. 

One  relative  aperture  distribution  considered  is  [1  - (p/a)2]2,  where  p 
is  the  radial -di stance  variable  in  the  aperture,  0 ^ p ^ a,  and  a = D/2. 

Figure  1 shows  this  distribution's  relative  power-density  contours  in  the  y-z 
plane  for  the  case  D > 30A,  where  -4D  ^ y ^ 4D  and  0 ^ z ^ D2/A.  Each  contour 
corresponds  to  increments  of  2.5  dB  below  the  peak  field  intensity;  the  darker 
contours  correspond  to  increments  of  5 dB  below  the  peak  field  intensity. 

With  a rotational ly  symmetric  aperture  distribution  the  contours  at 
positive  and  negative  y values  will  be  mirror  images.  Thus,  we  only  need  the 
range  0 ^ y ^ 4D.  Also,  using  this  range  doubles  the  accuracy  of  specifying 
the  graph  ordinate  y.  The  result  using  the  same  case  as  figure  1 is  shown  in 
figure  2.  Overlaying  the  contours  in  figure  2 are  dashed  lines  corresponding 
to  constant  y values,  with  incremental  spacing  of  ^ between  adjacent  lines. 

In  figure  3 the  relative  field  intensity  along  each  of  these  lines  is  plotted 
as  a set  of  parametric  curves.  The  uppermost  curve  in  figure  3 corresponds  to 

y I 0. 

The  preceding  figures  describe  the  case  D > 30A.  To  indicate  how  the 

results  change  for  smaller  ^ ratios,  in  figures  4,  5,  and  6 we  show  relative 

^ D 

y-z  plane  power-density  contours  for  T values  equal  to  24,  16,  and  10,  respec- 

A 

tively.  Relatively  little  difference  exists  between  figures  2 and  4,  whereas 
considerable  difference  exists  between  off-boresight  contours  in  going  from 
figure  2 to  figure  6.  For  further  comparison,  in  figures  7 and  8 we  show 
relative  field  intensity  curves,  along  equally-spaced  parallel  lines  analogous 
to  figure  3,  for  ^ values  of  16  and  10,  respectively. 

The  absolute  power  density  in  dBm  per  square  centimeter  at  any  point  in 
front  of  the  antenna,  for  one  watt  of  input  power  at  the  antenna  aperture,  is 
readily  arrived  at  by  combining  the  relative  dB  reading  obtained  from  the  pre- 
ceding figures  with  the  quantity,  38.57  - 20  x Log10D  [dBm/cm2],  where  D is 
in  centimeters. 
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We  next  present  results  for  a slightly  modified  relative  aperture  distri- 
bution, of  the  form  a + [1  - (^ ) 2 ] 2 , where  a is  a small  constant.  Thus,  in 

3 

figures  9 and  10  we  present  relative  power  density  contours  and  relative  field- 
intensity  curves  for  a = 0.1  (a  pedestal  20  dB  down),  while  in  figures  11  and 
12  we  present  results  for  a = 0.316  (a  pedestal  10  dB  down).  In  the  first 
case  the  absolute  power  density  in  dBm/cm2  for  one  watt  of  input  power  at  the 
antenna  aperture  is  arrived  at  by  combining  the  relative  dB  reading  obtained 
from  figure  9 or  10  with  38.17  - 20  x Log10D,  while  for  the  second  case  the 
absolute  power  density  is  obtained  by  combining  the  relative  dB  reading  from 
figure  11  or  12  with  38.14  - 20  x Log10D.  Again,  D is  in  centimeters. 

To  illustrate  the  use  of  these  graphs,  assume  that  the  antenna  is  a 11.3  m 

(37  ft)  diameter  parabolic  reflector  operating  at  6.5  GHz  with  56  kW  of  input 

power,  and  that  we  wish  to  estimate  the  power  density  823  m (2700  ft)  out  from 

the  antenna  face  and  15.2  m (50  ft)  off  axis.  Thus,  the  y distance  is  1.35D 

D2 

and  the  z distance  is  0.3^-.  From  figure  2,  we  see  that  at  these  coordinates 
there  is  a relative  power  density  of  -45  dB  relative  to  the  peak,  an  amount 
also  confirmed  by  figure  3.  The  correspondi ng  absolute  power  density  per  watt 
of  input  power  is  -67.5  dBm/cm2,  so  the  actual  power  density  is  about  0.01 
mW/cm2  at  this  location.  The  effect  here  of  a pedestal  in  the  relative 
aperture  distribution  is  destructive  interference  according  to  figures  9 and 
11,  so  the  above  estimate  may  be  considered  as  worst  case. 

In  conclusion  then,  we  have  presented  a nomograph  and  associated  algorithm 
for  obtaining  the  absolute  power  density  in  the  near  field  of  a circular  aper- 
ture with  a [1  - (£)2]2  relative  aperture  distribution,  along  with  results  for 
this  distribution  on  a 10  dB  pedestal  and  on  a 20  dB  pedestal.  It  should  be 
noted,  in  extrapolating  these  results  to  actual  antennas,  that  aperture  block- 
age and  struts  will  produce  additional  effects  that  have  not  been  accounted 
for  here. 
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Figure  1.  Relative  Power  Density  Contours  in  the  Y-Z  plane  for  a [1  - (^)  ] 
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aperture  distribution:  D > 30\;  -4D  to  +4D  range  on  Y. 
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axis  distance  divided  by  D 


Figure  2.  Relative  Power  Density  Contours  in  the  Y-Z  plane  for  a [1  - (^)2]2 

3 

aperture  distribution:  D > 30A. 
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Relative  Field  Strength  in  dB 


Z-axis  distance  divided  by  D2/A 


Figure  3.  Relative  Field  Intensity  along  equally  spaced  parallel  lines  for  a 
[1  “ (§)2]2  aperture  distribution:  D > 30A. 
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Figure  4.  Relative  Power  Density  Contours  in  the  Y-Z  plane  for  a [1  - (§)2]2 

a 

aperture  distribution:  D = 24A. 
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Figure  5.  Relative  Power  Density  Contours  in  the  Y-Z  plane  for  a [1  - (£)2]2 

a 

aperture  distribution:  D = 16A. 
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-axis  distance  divided  by  D 


Figure  6.  Relative  Power  Density  Contours  in  the  Y-Z  plane  for  a [1  - (^)2]2 
aperture  distribution:  D = 10A. 
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Relative  Field  Strength  in  dB 
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Figure  7.  Relative  Field  Intensity  along  equally  spaced  parallel  lines  for  a 
[1  - C^)2]2  aperture  distribution:  D = 16A. 
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Relative  Field  Strength  in  dB 


Z-axis  distance  divided  by  D2/A. 

Figure  8.  Relative  Field  Intensity  along  equally  spaced  parallel  lines  for  a 
[1  - (£)2]2  aperture  distribution:  D = 10A. 
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Figure  9. 


Relative  Power  Density  Contours  in  the  Y-Z  plane  for  a [1  - (^)2]2 

3 

aperture  distribution  on  a pedestal  20  dB  down:  D > 30A. 
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Relative  Field  Strength  in  dB 
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Figure  10.  Relative  Field  Intensity  along  equally  spaced  parallel  lines  for  a 
[1  - (^)2]2  aperture  distribution  on  a pedestal  20  dB  down:  D > 30A. 
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Figure  11.  Relative  Power  Density  Contours  in  the  Y-Z  plane  for  a [1  - (^)2]2 

3 

aperture  distribution  on  a pedestal  10  dB  down:  D > 30A. 
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Figure  12.  Relative  Field  Intensity  along  equally  spaced  parallel  lines  for  a 
[1  - (£)2]2  aperture  distribution  on  a pedestal  10  dB  down:  D > 30A. 
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III.  Computation  of  the  Electric  Field  Intensity  in  the  Near-Zone  of 
Aperture  Antennas  by  Numerical  Integration 
of  the  Plane-Wave  Spectrum 

Considerable  attention  has  been  given  to  the  problem  of  calculating  the 
far-field  pattern  of  an  antenna  from  measured  near-field  data  [1]  - [11].  Of 
these  techniques,  planar  near-field  scanning  (plane-wave  spectrum  deconvolu- 
tion) has  seen  the  most  development  [12]  - [16].  Here,  we  wish  to  consider 
the  inverse  problem  of  determining  near-zone  and  Fresnel  region  fields  from 
numerical  values  of  the  complex,  vector,  far-field  pattern.  That  is,  we  are 
interested  in  numerical  evaluation  of  the  plane-wave  spectrum  integral  repre- 
senting the  radiated  electric  field.  We  assume  that  numerical  plane-wave  spec- 
trum data  has  previously  been  generated  either  from  planar  near-field  measure- 
ments or  from  theoretical  far-field  patterns. 

An  earlier  algorithm  was  developed  by  Yaghjian  [17]  for  carrying  out  such 
a numerical  integration.  In  Yaghjian1 s analysis,  the  integration  domain  was 
restricted  in  order  to  achieve  computational  efficiency;  this  also  served  to 
constrain  the  oscillatory  nature  of  the  integrand  and  so  eliminate  a cata- 
strophic numerical  instability  which  otherwise  occurs  at  large  z-axis  dis- 
tances with  discrete  data  points.  Use  of  the  fast  Fourier  transform  (FFT) 
provides  efficient  numerical  integration  over  a large  transverse  output  range. 

Extending  this  earlier  work  to  measured  data  with  fixed  data-point  spacing 

D2 

and  to  z-axis  distances  greater  than  0.5  — , where  D is  the  antenna  diameter 
and  A is  the  wavelength,  we  found  that  numerical  instabilities  appeared  at  the 
extremities  of  the  integration  range.  Attempts  to  reduce  these  instabilities 
by  uniformly  decreasing  the  already  restricted  integration  domain  drove  the 
calculated  field  to  zero  at  the  extremities,  an  effect  that  had  also  been  found 
by  Yaghjian  [17].  In  order  to  compensate  for  these  computational  difficulties, 
a new  limiting  criterion  was  developed. 

In  Yaghjian1 s analysis,  a geometrical  (sheaf  of  angles)  criterion  was 
used  to  limit  the  integration  domain.  In  this  work,  we  develop  an  improved 
formula  for  limiting  the  integration  domain,  using  a sampling-theorem  [18] 
criterion,  producing  a numerically  more-stable  result  over  a greater 
cartesian-coordi nate  output  range.  Also,  we  establish  upper  bounds  on  the 
z-axis  distance  from  the  antenna  and  on  the  data-point  spacing  interval  (A)  of 
the  plane-wave  spectrum  in  order  to  ensure  reliable  computations. 

Evaluating  the  plane-wave  spectrum  integral,  then,  we  compute  near-zone 
and  Fresnel-region  fields  at  a sequence  of  fixed  z-axis  distances  as  a 
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function  of  the  lateral  coordinates  x and  y.  We  reduce  computational  effort 
by  alternately  setting  one  lateral  coordinate  equal  to  zero.  By  computing 
along  just  the  x and  y axes  we  obtain  sufficient  data  to  generate  near-zone 
field-intensity  contours  in  the  x-z  and  y-z  planes.  The  maximum  lateral  coor- 
dinate range  is  Rq  = However,  we  empirically  determined  that  an  effective 
lateral  coordinate  range  about  half  this  size  is  needed  to  avoid  numerical 
instabilities.  Thus,  this  effective  lateral  coordinate  range  and  the  upper 
bound  on  the  z-axis  distance  together  define  a truncated  cylindrical  region 
bounding  the  near-field  computation  domain. 

Yaghjian  [17]  describes  how  a coordinate-system  rotation  of  the  plane- 
wave  spectrum  data  allows  the  cylindrical  region's  axis  to  subtend  an  arbi- 
trary angle  with  respect  to  the  main  beam  of  the  antenna;  however,  this  also 
results  in  havi:.g  to  interpolate  when  using  discrete  plane-wave  spectrum  data. 
Here,  by  aligning  the  cylindrical  region's  axis  with  the  perpendicular  to  the 
near-field  measurement  plane,  interpolation  becomes  unnecessary  and  so  we 
avoid  interpolation  error. 

To  be  definite,  we  describe  the  case  where  the  plane-wave  spectrum  data 
are  obtained  from  planar  near-field  measurements.  The  antenna  measurement 
setup  is  depicted  in  figure  13,  showing  the  antenna  under  test,  the  near-field 
measurement  plane  with  the  probe  antenna  just  in  front  of  the  test  antenna, 
and  a hypothetical  y-z  plane  cut  through  the  cylindrical  region  of  interest. 
Graphical  output  is  represented  in  the  figure  by  field-strength  contours  in 
the  y-z  plane.  At  equally  spaced  intervals  along  the  z-axis,  vertical  dashed 
lines  are  drawn  to  represent  the  series  of  y-axis  cuts  along  which  the  near 
field  is  computed  by  the  planar  near-field  program.  For  smooth  contours,  of 
course,  many  more  plane-cut  computation  points  would  be  required  along  the 
z-axis  than  are  indicated  by  the  figure. 

As  background,  and  to  introduce  notation,  we  give  an  expression  for  the 
y-component  of  a transmitting  antenna's  plane-wave  spectrum  in  terms  of  mea- 
sured probe  data  taken  over  a planar  surface  in  front  of  the  antenna.  It  is 
assumed  that  multiple  reflections  between  the  transmitting  antenna  and  the 
probe  are  negligible,  that  both  antennas  are  polarization  matched  in  the  same 
(y)  direction,  that  the  input  reflection  coefficients  at  the  waveguide  leads 
to  the  antennas  are  negligible,  and  that  the  probe  is  far  enough  removed  from 
the  transmitting  antenna  so  that  the  significant  probe  receiving  pattern  is 
independent  of  spatial  frequency  variations,  i.e.,  the  incident  radiation  at 
the  probe  has  a relatively  narrow  angle  of  arrival.  Let  B (P)  be  the  signal 
received  by  the  probe  at  position  P on  the  planar  scan  surface  normalized  to 
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intensity  Contours  are  Plotted. 


unity  at  position  P1  (here,  P = xa  + ya  , where  a , a , a are  unit  cartesian 

x y x y z 

vectors),  and  let  Aq  be  the  insertion  loss  between  the  transmitting  antenna 
and  the  probe  when  it  is  located  at  position  P1 . Then  the  plane-wave  spec- 
trum's y-component  can  be  expressed  [12,13]  as 


b(K) 


1 f 2Po 

\ V nAYG(o) 
^ o o r 


B (P)e"1-*-  dP, 
o — — 


(1) 


where  Pq  is  the  input  power  to  the  antenna,  Yq  is  the  characteristic  admit- 
tance of  free  space,  G^Co)  is  the  gain  of  the  probe,  d is  the  z-axis  distance 
from  the  planar  scan  surface  to  the  transmitting  antenna,  dP  = dxdy,  A.  is  the 
wavelength,  K = k a + k a , and  y = , jk 2 - K2  where  k = 2n/\  and  K2  = K-K. 

The  function  Bq ( P ) is  virtually  representable  (i.e.,  b ( K)  approaches  zero 
for  K sufficiently  greater  than  k,  depending  upon  d)  as  the  Fourier  transform 
of  a bandlimited  function.  Consequently,  the  sampling  theorem  applies,  so  the 
double  integral  (1)  is  completely  determined  by  the  sampling  of  data  at  dis- 
crete equal-spaced  lattice  points  over  the  scan  plane  [12].  If  we  assign  band 

limits  ±k  and  ±k  to  k and  k , respectively,  then  (1)  becomes  [ibid] 
xm  ymx  y 
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To  be  mathematically  exact,  N2  and  N2  should  be 


infinite;  in  practice,  Nx  and  N2  are  limited  by  the  size  of  the  scan  area.  As 
given  above,  the  FFT  would  limit  to  2Ni  x 2N2  the  number  of  points  at  which 
values  of  b ( K)  would  be  calculated.  However,  we  increase  the  output  spectral 
density  (number  of  data  points)  by  zero  filling  around  the  scan  area  boundary. 
Thus,  setting  BQ(Prs)  = 0 for  |r|  > Nx  or  for  |s|  > N2 , we  have 
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where  - ^ m ^ - 1,  * ^ n ^ • 1.  Here,  Kmn  = mAxax  + nAyay 

2 k 2 k 

y = Vk2  “ K2  , where  A = Klx-  , A = -jr—  and  N > 2Nj , N > 2N2. 
Jmn  N mn’  xN’yN  x y 


and 
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We  now  return  to  the  development  of  the  near-field  computation  algorithm. 

The  electric  field  can  be  obtained  at  any  point  r = R + zaz,  where 

R = xa  + ya  , from  the  plane-wave  spectrum  integration 
x y 


E(r)  = ~-f f b(K)e'-"-dK 

-OO 


(3) 


where  k = K + ya  and  dK  = dk  dk  . Here  we  assume  that  the  vector  spectral- 
density  function  b(K)  is  known.  If  both  the  x and  y components  of  b(K)  have 
been  determined,  then  the  z component  is  obtained  from  b^ ( K)  = - ^ K*b(K).  In 
the  far-field  of  the  antenna  (3)  becomes 


, i kr 

E(r)  ~ -ik  cos0  b(£  R) 


where  cos0  = z/r  and  r = y[r*r  . 

In  the  Fresnel  region  evanescent  modes  will  not  be  significant,  so  the 
integration  limits  in  (3)  change  from  ±°°  to  ±k.  The  minimum  applicable  z dis- 
tance is  at  least  a few  wavelengths  to  ensure  that  evanescent  modes  have  at- 
tenuated. In  our  technique,  the  integration  is  carried  out  using  the  FFT  for 
a single  fixed  value  of  z at  a time.  However,  stability  problems  can  become 
critical  in  carrying  out  a numerical  integration.  In  particular,  when  z is 
moderately  large,  the  factor  e1^2  in  the  integrand  can  oscillate  rapidly 
between  plus  and  minus  one  while  the  rest  of  the  integrand  changes  very 
little.  As  a result,  these  terms  cancel  out  analytically,  but  with 
discrete  data  on  a computer  they  can  add  numerically  to  produce  a significant 
integration  error.  This  can  occur  when  z/D  is  as  small  as  2 or  3,  where  D is 
the  antenna  diameter.  One  can  compensate  for  this  by  limiting  the  domain  of 

integration  even  further,  say  to  -k  ^ k ^ k ,-k  ^ k ^ k , where  k < 

3 J xo  x xo’  yo  y yo’  xo 

k,  kyQ  < k.  That  is,  we  have  to  evaluate  the  near-field  expression 


E(r) 
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(4) 
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We  obtain  integration-domain  limits  for  (4)  by  using  sampling  theorem 
criteria;  namely,  the  maximum  change  in  yz  between  two  adjacent  sample  points 
shall  be  less  than  7T.  After  some  ado,  this  leads  to  the  result 


*S<o  * *Vo  _ Ro  (5) 

k ^/fs+z2 

where  R s n/ A E n/A  . Here,  B = 1 + sin  2d>  where  d>  = tan  "k  /k  . The 
o y x ’ M Yo  Yo  yo  xo 

derivation  of  (5)  is  presented  in  the  appendix,  where  we  establish  that  the 

R2 

maximum  z distance  is  limited  by  z < 

Using  the  FFT  to  evaluate  (4),  it  is  expedient  to  integrate  over  the 

extended  domain  |k|^k  ,|k|^k  , where  k = k ^ k,  and  utilize  a 

x xm  ’ y ym  ’ xm  ym  ’ 

window  function  to  limit  the  effective  integration  domain  to  k < k , 

a xo  xm’ 

kyo  < ^ym'  ^hus,  for  a fixed  value  of  z,  we  can  approximate  (4)  by 


Here,  E.  0 is  the  Fresnel-region  electric  field  in  the  x-y  plane  corresponding 

to  a fixed  value  of  z with  x = j6  , y = £6  , where  6 = zt/k  and  6 = n/k 

J x’  J y’  x xm  y ym 

correspond  to  measurement-plane  data-point  spacing  increments.  Also,  F(kx,ky) 

= b(K)u(K)e1^z,  where  u(K)  is  a window  function  for  limiting  the  oscillatory 

behavior  of  e1^2.  The  extended  integration  domain  limits  are  given  by  kxm  = 

1/2  N A , k = 1/2  N A . In  (6),  -N  /2^j<N  /2,  -N  /2^£<N  /2,  so  the  maximum 
xx’  ym  yy  v ’ x J x ’ y y ’ 

lateral  coordinate  range  is  |x|^l/2  N 6 = n/A  , |y|^l/2  N 6 = n/A  . How- 

xx  x y y y 

ever,  good  numerical  stability  of  (6)  at  moderately  large  z values  limits  the 
effective  lateral  coordinate  range  to  about  | x | < 1/4  Nx<5x,  | y | < 1/4  N 6 . 

Equation  (5)  was  compared  against  an  alternate  expression  obtained  by 
Yaghjian  [17].  Although  (5)  yielded  the  best  results,  it  lacked  the  latter's 
ease  of  implementation.  However,  with  a hybrid  formula,  we  obtained  an  ex- 
pression that  yielded  good  results  along  with  ease  of  computation.  According- 
ly, the  window  function  u(K)  is  set  equal  to  zero  under  the  condition 
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Ikxl  ^ N > Ro  + D 
k 

y o 

The  integration  cutoff  imposed  by  u(K)  can  lead  to  ringing  at  the  extremities 
of  the  computed  near  field  when  the  amplitude  of  b(K)  at  the  cutoff  point  is 
large.  This,  in  fact,  is  the  reason  for  limiting  the  effective  range  on  x and 
y to  about  half  the  normal  FFT  output  range.  Note  that  this  range  limitation 
quadruples  the  number  of  data  points  needed  for  the  plane-wave-spectrum 
i ntegration. 

We  now  write  (6)  in  a form  more  suitable  for  FFT  processing.  The  result 
i s 


N -1  N -1 


E j ‘ * 1 ' 6 J ] m=o  Jo  {- (m  ' - >y3e^B)")  exp  + £L)}  , 

c X y 


where  the  ranges  0^j^Nx~l,  0S£^Ny-l  apply.  As  a result,  our  formulation  agrees 
with  standard  FFT  practice.  It  is  assumed  that  Nx/2  and  N^/2  are  both  even 
i ntegers. 

In  case  we  are  interested  in  the  field  along  the  y-axis  where  x = 0,  the 
preceeding  expression  becomes 
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With  this  formulation,  we  simply  sum  along  the  k coordinate  followed  by  a 
one-dimensional  FFT  on  the  coordinate.  A similar  "one-dimensional  col- 
lapse" formulation  is  used  to  compute  the  near  field  along  the  x-axis. 

Near-field  computations  using  this  algorithm  were  tested  against  both 
theoretical  calculations  and  experimental  data.  Thus,  in  figures  14  and  15  we 
compare  results  produced  by  the  planar  near-field  program  against  an  exact 
theoretical  expression  obtained  by  Rudduck  and  Chen  [19]  for  the  field  along 
the  z-axis  of  a uniform  circular  aperture.  The  absolute  power  densities  are 
obtained  by  each  method  and  compared  in  the  figures  for  both  a 10A  diameter 
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circular  aperture  and  a 20A  diameter  circular  aperture  at  an  operating 
frequency  of  4 GHz.  Computed  results  are  plotted  as  solid  lines  while 
theoretical  results  are  plotted  as  dashed  lines.  The  on-axis  electric-field 
strength  for  a uniformly  excited  circular  aperture  is  given  by 


where  a = D/2  is  the  aperture  radius.  The  close  correspondence  between  com- 
puted and  theoretical  results  in  these  figures  is  extremely  gratifying. 

Next,  in  figures  16  through  21,  we  show  comparisons  between  computed 
near-zone  results  and  actual  measured  data  taken  at  corresponding  z-axis 
distances.  The  antenna  characterized  in  figures  16  through  18  (antenna  #1) 
has  D/A  ratio  of  16.2,  while  the  antenna  characterized  in  figures  19  thru  21 
(antenna  #2)  has  a D/A  ratio  of  77.6.  For  each  antenna,  we  compare  measured 
data  against  predicted  near  fields  at  the  measurement  plane  and  at  z-axis 
distances  of  100  cm  and  300  cm  beyond  that.  Both  antennas  are  linearly 
polarized  in  the  y-di recti  on;  also,  the  reconstructed  near-fields  and  the 
comparison  measurements  are  taken  parallel  to  the  y-axis.  As  seen  in  the 
figures,  there  are  fewer  reconstructed  near-field  points  than  measurements  at 
the  shortest  distance,  due  to  using  a smaller  scan  area  for  generating  the 
plane-wave  spectrum  than  for  taking  comparison  measurements.  Measured  results 
are  plotted  as  dashed  lines  while  the  computed  results  are  plotted  as  solid 
lines.  Both  measured  values  and  computed  results  are  normalized  with  respect 
to  the  corresponding  peak  field  strength  in  the  measurement  plane.  The 
agreement  shown  between  measured  and  computed  results  is  remarkable  in  view  of 
the  very  complex  pattern  structure  of  these  near-zone  fields. 

Finally,  in  figures  22  and  23  we  show  computed  near-zone  power-density 
contours  in  the  y-z  plane  for  these  same  two  antennas,  while  in  figure  24  we 
show  theoretical  near-zone  power-density  contours  for  a 10A  diameter 
uniformly  excited  circular  aperture  at  4 GHz.  These  contours  give  absolute 
power  density  in  dBm/cm2  corresponding  to  one  Watt  of  input  power  at  the 
antenna  aperture. 

In  conclusion,  we  have  presented  an  improved  algorithm  for  calculating 
near-zone  fields  in  the  vicinity  of  radiating  aperture  or  dish  antennas  and  we 
have  demonstrated  that  our  algorithm  is  capable  of  extremely  high  accuracy. 

We  have  also  presented  absolute  power-density  contours  for  some  of  the 
antennas  studied.  Such  absolute  power-density  contours  are  a convenient  and 
informative  method  for  exhibiting  near-zone  computation  results. 
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Appendix 


Here,  we  wish  to  obtain  integration  limits  kxQ,  k^o  to  bound  the  oscil- 
latory behaviour  of  e]^z  in  the  integrand  of  (4).  We  do  this  by  making  use  of 
the  sampling  theorem,  which  requires  that  there  be  at  least  two  sample  points 
per  cycle  at  the  greatest  resolvable  "z"  spatial  frequency.  That  is,  we  will 
only  retain  those  z-component  spatial  frequencies  within  the  integrand  satis- 
fying this  criterion,  so  the  maximum  allowable  incremental  change  in  y will  be 
71/z.  Since  y = Vk2  - K2 , where  K2  = k£  + k2,  the  applicable  criterion  for 
determining  the  integration  domain  limits  is 


L = V 1 - (U*  - b)  - V 1 - U2  , 


(8) 


where  U = -r—  = V k2  + k2  /k,  and  where  b represents  a small  incremental 
k ^ xo  yo  r 

change  in  U2  corresponding  to  a change  from  the  sample  point  kxQ,  k^Q  to  the 

point  k - A , k - A . For  consistency  we  require  that  k + k > 
r xo  x’  yo  y J ^ xo  yo 

A + A . 


x y 


The  assumption  that  evanescent  modes  are  negligible  in  (3)  implies  that 

j U2 

U < 1.  Accordingly,  we  can  apply  the  binomial  expansion,  1 - U2  - 1 - r - 

U4  l 

g-  - to  (8).  After  cancelling  terms  and  rearranging  series  expansions  we 

obtai n 


7i  _ b /-j  , 1 2 U8  5U10 

kz  8 y1  1 - U2  1 - U2/2  32  64 


(9) 


Neglecting  terms  of  order  b2  or  smaller,  we  can  express  b as  the  ratio  of  a 
fourth  order  polynomial  to  a power  series  in  U2.  Carrying  out  the  indicated 
power-series  division  we  obtain 
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(10) 
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Now  let  A e A = A . Then 

x y 

k2b  = K2  - (k  - A)2  - (k  - A)2  = 2A  (k  + k - A) 
o xo  yo  xo  yo 


or  kb  = 2A(aU  - ==) , where  a = cos  <j)  + sin  4>  ; 4>  = tan 


-i  i\ 

■l  yo 
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We  have  (see  text  following  eq.  2)  A = n/R  , where  R = 6 is  the  maximum 

o o2yy 

lateral-coordinate  range.  Consequently,  b2  in  (9)  is  negligible  compared  to  b 
provided  that  Rq  >>  A.  Now  we  can  express  (10)  as 


U2  + 2 jp  U 
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Az  _ 0 
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(11) 


We  solve  (11)  as  a simple  quadratic  equation 

A 


. U4  U6  5U8 

where  o = -j  + -g  + -p  + 

in  U.  Then,  for  z > Rq  and  neglecting  terms  smaller  than  , we  expand  the 

resulting  radical  as  a power  series  in  ^ including  the  small  correction  term 

ct  in  the  expansion.  Iteratively  substituting  into  the  expression  for  a we 

/ R \13 

obtain  a corrected  series  expansion  for  U,  which  to  an  accuracy  of  — ) is 


az 


equivalent  to 


U = 
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f + 1 + h* 

R 1 R2 
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(12) 


Equation  (12) 's  greatest  source  of  discrepancy  comes  from  neglecting  the  b2 

AR 

terms  in  (9),  so  for  z large  the  discrepancy  in  U is  on  the  order  of  r— 2 

y / P / x \2 

Thus,  for  p-  >(_o^2  ''  1 

o V A 

see  that  our  consistancy  condition,  k + k > 2A,  is  satisfied  provided 

J ’ xo  yo  ’ K 


, the  discrepancy  is  less  than( p- ) . From  (12)  we  readily 


„ 2R2 
z < o 

A 


Finally,  neglecting  terms  of  order  A/R  , (12)  becomes 


U = 


r' 
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Absolute  Power  Density  in  dB  below  one  mi  1 1 i -Watt/cm 
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Figure  14.  Comparison  of  results,  along  the  z-axis,  as  produced  by  the  planar 
program  (solid  line)  and  as  calculated  using  an  exact  theoretical 
expression  (dashed  line)  for  a 1 0 — A diameter  uniformly  excited 
circular  aperture  at  4 GHz  showing  power  density  in  dB  below  one 
milliwatt/cm2,  assuming  one  watt  input  power  to  the  aperature. 
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Absolute  Power  Density  in  dB  below  one  mi  1 1 i -Watt/cm 
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Figure  15.  Comparison  of  results,  along  the  z-axis,  as  produced  by  the  planar 
program  (solid  line)  and  as  calculated  using  an  exact  theoretical 
expression  (dashed  line)  for  a 20-A  diameter  uniformly  excited 
circular  aperture  at  4 GHz  showing  power  density  in  dB  below  one 
milliwatt/cm2,  assuming  one  watt  input  power  to  the  aperature. 


30 


Relative  Field  Strength  In  dB 


Y-Axis  Distance  divided  by  D 


Figure  16.  Comparison  of  Measured  (dashed  line)  and  Computed  (solid  line) 
Relative  E-Plane-Cut  Field  Intensity  Curves  for  Antenna  #1  at 
z = 65.3  cm. 
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Figure  17.  Comparison  of  Measured  (dashed  line)  and  Computed  (solid  line) 
Relative  E-Plane-Cut  Field  Intensity  Curves  for  Antenna  #1  at 
z = 165.3  cm. 
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Figure  18.  Comparison  of  Measured  (dashed  line)  and  Computed  (solid  line) 
Relative  E-Plane-Cut  Field  Intensity  Curves  for  Antenna  #1  at 
z = 365.3  cm. 
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Figure  19.  Comparison  of  Measured  (dashed  line)  and  Computed  (solid  line) 
Relative  E-Plane-Cut  Field  Intensity  Curves  for  Antenna  #2  at 
z = 42  cm. 
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Figure  20.  Comparison  of  Measured  (dashed  line)  and  Computed  (solid  line) 
Relative  E-Plane-Cut  Field  Intensity  Curves  for  Antenna  #2  at 
z = 142  cm. 
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Figure  21.  Comparison  of  Measured  (dashed  line)  and  Computed  (solid  line) 
Relative  E-Plane-Cut  Field  Intensity  Curves  for  Antenna  #2  at 
z = 342  cm. 
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Figure  22.  Near-Zone  Power  Density  Contours  for  Antenna  #1  in  dB  below  one 
mil  1 i-Watt/crrr  for  one  Watt  input  power  to  the  antenna  (2.5  dB 
separation  between  adjacent  contours). 
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Figure  23.  Near-Zone  Power  Density  Contours  for  Antenna  #2  in  dB  below  one 
mi  1 1 i-Watt/cnr  for  one  Watt  input  power  to  the  antenna  (2.5  dB 
separation  between  adjacent  contours). 
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Figure  24.  Near-Zone  Power  Density  Contours  for  a 10-A  Diameter  Uniformly 

Excited  Circular  Aperture  in  dB  below  one  mi  11 i -Watt/cm  at  4 GHz 
for  one  Watt  input  power  at  the  aperture  (2.5  dB  separation  betweer 
adjacent  contours. ) 
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